1 


NASA 

Technical 

Paper 

3069 


March 1991 


NASA-TP-3069 WlOOO^HH O 


Structural Factoring 
Approach for Analyzing 
Stochastic Networks 


Kelly J. Hayhurst 
and Douglas R. Shier 



n nr* r* n 


Lewsy.ih c uu 



•SiEV RESEARCH CENTER 
LIBRARY flASA 
HAMPTON, VIRGINIA 







NASA 

Technical 

Paper 

3069 


1991 


Structural Factoring 
Approach for Analyzing 
Stochastic Networks 


Kelly J. Hayhurst 
Langley Research Center 
Hampton , Virginia 

Douglas R. Shier 

The College of William and Mary 

Williamsburg, Virginia 


NASA 

National Aeronautics and 
Space Administration 

Office of Management 

Scientific and Technical 
Information Division 



IV 


Contents 

Symbols 

Abstract 

Introduction 

Notation 

Structural Reduction Techniques 

Fundamental Reductions 

Conditioned Reductions 

Conditioned Reductions and Shortest Path Calculations 

Total Factoring 

Determining the Shortest Path Length Distribution 

Applications of the Structural Factoring Approach 

Critical Path Analysis 

Two-Terminal Reliability 

Comparison Between Structural Factoring and Complete Enumeration 

Conclusions 

References 


1 

1 

2 

2 

3 

5 

10 

13 

13 

15 

15 

16 
18 

18 

19 


iii 



Symbols 

A 

a (tJ) 

c j 

/(<) 

fi(j) 

G 

(iJ) 

indegree(i) 

771 

N 

71 

outdegree(z) 

PQ 

P (hj) 

Sj 

S 

t 

7 

u(i) 

1 (U)\ 


set of arcs in a directed graph 
n-tuple of lengths for arc (z, j) 

j th length in an n-tuple 
factoring arc of node i 
zth factoring arc of node j 
directed graph 

arc with origin node i and destination node j 

number of arcs entering node i 

constant length for arc (z, j) 

particular length for arc (z, j) 

number of arcs in the network 

set of nodes in a directed graph 

number of nodes in the network 

number of arcs leaving node z 

probability function 

probability distribution for the lengths of arc (z, j) 

subnetwork j 

source node 

terminal node 

number of factoring arcs 

number of subnetworks generated by factoring on node z 
number of lengths for arc (i,j) 


iv 



Abstract 

In this report, the problem of finding the distribu- 
tion of the shortest path length through a stochas- 
tic network is investigated. A general algorithm is 
developed for determining the exact distribution of 
the shortest path length. The algorithm is based 
on the concept of conditional factoring, in which a 
directed, stochastic network is decomposed into an 
equivalent set of smaller, generally less complex sub- 
networks. Several network constructs are identified 
and exploited so that the computational effort re- 
quired to solve a network problem is significantly less 
than that required by complete enumeration. This 
algorithm can be applied to two important classes 
of stochastic path problems: determining the crit- 
ical path distribution for acyclic networks and the 
exact two-terminal reliability for probabilistic net- 
works. Computational experience with the algorithm 
has been encouraging and has allowed the exact so- 
lution of networks previously analyzed only by ap- 
proximation techniques. 

Introduction 

In general, network analysis involves the study of 
systems that can be modeled in terms of nodes and 
arcs connecting certain pairs of nodes. Quantitative 
information, such as length, reliability, or time to 
completion, may also be associated with the arcs. For 
example, figure 1 shows a simple six-node network 
representing possible shipping routes from a factory 
(node A) to a warehouse (node F). The arcs of 
this network represent highway routes with their 
associated mileage given as the arc lengths. 


Figure 1. Elementary shipping network. 

A commonly encountered problem in such a net- 
work is finding the shortest path from the factory 
to the warehouse. More generally, shortest path cal- 
culations are found to be valuable in analyzing the 
behavior of various large-scale distribution networks, 
such as complex telecommunications processes, dis- 
tributed computer architectures, and even lifeline 
systems subject to seismic risk (ref. 1). The im- 
portance of these applications requires that efficient 
techniques and tools be developed to aid in the net- 
work analysis process. 


Although techniques for deterministic network 
analysis are widespread and applicable to large-scale 
networks (ref. 2), this is not the case for the analysis 
of stochastic networks, where the behavior of the net- 
work components is governed by some random pro- 
cess. The most straightforward approach for investi- 
gating the stochastic shortest path problem involves 
enumeration of all possible states of the network. In 
general, if the network contains m arcs and each arc 
can take on k values, there would be k m states of the 
network that must be analyzed. This approach is ob- 
viously computationally infeasible even for relatively 
small networks. Thus, current analysis methods have 
been largely confined to approximation, simulation, 
and bounding techniques (refs. 3 to 7). 

Since the utility of a network model critically de- 
pends on the suitability of its assumptions, a general 
model that does not impose restrictions on the dis- 
tributional form of the network components is desir- 
able. This report considers a broad class of stochastic 
networks with the following characteristics: (1) arc 
values (length, duration, cost, etc.) are discrete ran- 
dom variables, (2) these random variables are sta- 
tistically independent, and (3) nodes do not fail. 
These assumptions greatly increase the tractability 
of the problem while preserving the modePs realism 
by incorporating the random nature of the problem 
(ref. 8, p. 455). Statistical independence is commonly 
assumed because it reduces the computational re- 
quirements associated with most solution approaches. 
Additionally, there are no assumptions, such as sym- 
metry of the nodes or arcs, imposed on the architec- 
ture of the networks. 

This report considers the specific problems of 
finding the distribution of the length of the shortest 
or longest (i.e., critical) path through such a stochas- 
tic network. The arcs of the network assume random 
lengths, which could represent, for example, the du- 
ration in traversing a communication link or the cost 
of completing a given function. In real-world set- 
tings, these values are most realistically viewed as 
random variables, rather than as fixed determinis- 
tic parameters. Thus, the state of the network de- 
pends on the state assumed by each arc. It follows 
that the shortest (or longest) path through the net- 
work is a function of the random arc lengths; hence, 
path length can be characterized by a probability 
distribution. 

The objective of this report is to present a gen- 
eral algorithm for determining the exact distribu- 
tion of the shortest path length in a directed, sto- 
chastic network. This algorithm has applications to 
two important classes of stochastic path problems: 
determining the longest, or critical, path distribu- 
tion for acyclic networks and the exact two-terminal 




reliability for probabilistic networks. The structure 
of the paper is as follows. In the next section, the 
necessary notation and basic terminology are defined. 
The structural decomposition technique is described 
in the third section, and it is applied to solving the 
critical path and two-terminal reliability problems 
in the fourth section. Several network examples are 
given in this section. A summary of the findings in- 
cluding the limitations of the approach is presented 
in the last section. 

Notation 

To facilitate the discussion of probabilistic net- 
works, some basic terminology and notation are first 
introduced. A probabilistic network is modeled us- 
ing a directed graph G = (N, A), where N is a set 
of nodes, representing, for example, warehouses or 
communication centers in the network, and A is a 
set of arcs, representing traffic routes or communica- 
tion buses, connecting certain pairs of nodes. An arc 
(i, j), where i and j are elements of AT, is defined to 
be a directed link from the origin node i to the desti- 
nation node j. If more than one arc connects a pair of 
nodes, the arcs are denoted with different numbered 
superscripts. For example, the two arcs spanning 
nodes D and E in figure 2 are denoted (D^) 1 and 
(D,E) 2 . Since we are concerned only with the short- 
est (or longest) path between two given nodes, we 
consider only graphs with one source node denoted 
by s, which has only outgoing arcs, and one terminal 
node denoted by t , which has only entering arcs. Fig- 
ure 2 shows a directed graph with source node A and 
terminal node E. The indegree of node i, indegree(z), 
is the number of arcs entering node z, and the outde- 
gree of i, outdegree(i), is the number of arcs leaving i. 
For each of the networks considered, indegree(s) = 0 
and outdegree(tf) = 0. In figure 2, indegree(C) = 2 
and outdegree(C) = 1. 



A path in the graph from node i to node j is 
defined as an ordered sequence of arcs connecting 
the two nodes. For example, one possible path 
from node A to node E in figure 2 is through arcs 
(A,B) and (B,E). A cycle is a special type of path 


connecting a node to itself. In figure 2, arcs (C,B), 
(B,D), and (D,C) form a cycle. In general, the 
networks treated here are allowed to contain cycles 
except in the special case of longest path calculations. 
In calculating the length of the longest path through 
a network, the underlying graph must be acyclic since 
cycles in a graph would lead to a longest path of 
infinite length. 

To incorporate information about the random be- 
havior of the arcs, each arc is assigned a finite n - tuple 
of nonnegative integer values, indicating, for exam- 
ple, lengths or durations. The n-tuple of arc lengths 
for arc (i,j) is denoted a^j), and the number of 
lengths associated with arc (i,j) is denoted |(i,j)|. 
For each arc (i, j), there is also a corresponding dis- 
crete probability distribution for the arc lengths and 
this is denoted Pfaj)- A particular length assumed 
by arc (i, j) is denoted L(i, j). Suppose the length 
of arc (B,D) in figure 2 is uniformly distributed on 
the interval [6,9]. Then the discrete n-tuple of arc 
lengths is represented by |(B,D)| = 4 and G(b,D) — 
< 6, 7, 8, 9 >. The corresponding probability distri- 
bution is given by P(b,D) = < 0.25, 0.25, 0.25, 0.25 >, 
and P[L(b,d) = 8] = 0.25. Note that using dis- 
crete probability distributions significantly increases 
the tractability of the problem in comparison with 
continuous distributions. Methods for discretizing a 
continuous distribution, as shown in Dodin (ref. 3), 
are available when the arc behavior is described by a 
continuous distribution. 

Structural Reduction Techniques 

In this section, a graph-theoretic procedure for 
determining the distribution of the shortest path 
length through a stochastic network is described. 
The objective of the approach is to apply certain 
conditioned reductions to a given network until the 
network is reduced to an equivalent (with respect to 
finding the distribution of the shortest path length) 
network having only two nodes. The arc connecting 
the two nodes in this reduced network provides the 
distribution of the shortest path length in the original 
network. 

In this report, two classes of structural reductions 
are applied to stochastic networks. The first class of 
reductions represent fundamental simplifications to 
the topology of the network including the series and 
parallel reductions. The second class of reductions 
are based on the concept of local “factoring.” The 
basic difference between the two classes is that the 
fundamental reductions yield one smaller equivalent 
network. Local factoring, on the other hand, gen- 
erates an equivalent set of subnetworks to be solved. 
Local factoring uses the divide and conquer ideology: 
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the original problem that is difficult to solve is de- 
composed into successively smaller problems that are 
easier to solve. Each type of reduction is discussed 
more fully in the following subsections. 


a (A C) “ <2,3,4, 5, 6, 7, 8, 9 > 

® : ►© 

P(A G) = <05,. 1 ,.15,.2,.2,.15,.1,.05> 


Fundamental Reductions 


Figure 4. Reduced series construct. 


The fundamental reductions can be applied to 
three basic constructs present in a network: series, 
parallel, and figure-eight. The term construct is 
used throughout this discussion to describe specific 
configurations of nodes and arcs that exist within a 
network architecture. Each fundamental reduction 
simplifies the given network graph by decreasing the 
number of nodes or the number of arcs (or both). 

One of the most easily recognizable constructs 
within a directed graph is the series construct. A 
series construct exists when some node B in the 
graph is the intermediate node between exactly two 
other nodes; that is, when indegree(B) = 1 and 
outdegree(B) = 1. Figure 3 gives an example of 
a series construct in which the discrete distribution 
function on arc (A,B) is uniform on [0,3] and the 
discrete distribution function on arc (B,C) is uniform 
on [2,6]. 


a (A,B) “ <0 ’ 1 ’ 2>3> a (B,C) “ < 2 >3,4, 5, 6> 

® kd *>© 

P(A,B) = <-25,.25,.25,.25> P(b,c) = <-2,.2,.2,.2,.2> 

Figure 3. Series construct. 


Since any path through node B must include arcs 
(A,B) and (B,C), these two arcs can be replaced 
by a single arc between nodes A and C. The set 
of lengths for the replacement arc (A,C) is the set 
of all possible additive combinations of arc lengths 
chosen from «(a,B) an ^ a (B,C)* F° r our example, 
L( A, C) = 3 occurs with a positive probability since 
the combinations L( A, B) = 0 and L( B, C) = 3 and 
L(A,B) = 1 and L(B,C) = 2 are possible states of 
that construct. Moreover, because of the statistical 
independence of arc lengths associated with distinct 
arcs, 


P[L(A, C) = 3] = P[L(A, B) = 0]P[L(B, C) = 3] 

+ P[L(A,B) = 1]P[L(B,C) = 2] 
= (0.25)(0.2) + (0.25)(0.2) = 0.1 

More generally, those lengths in the set G(a,C) that 
occur in several ways are denoted only once in G(a,C) 
and the corresponding probability is the sum of the 
contributing probabilities. The reduced construct, 
which corresponds to the discrete convolution of the 
arc distributions on (A,B) and (B,C), is shown in 
figure 4. 


Two or more arcs that connect the same pair of 
nodes constitute a parallel construct. Figure 5 shows 
a parallel construct where for arc (AjB) 1 , G(ab) 1 = 
< 2,3,5 > and P(A,B)i = < 0.25, 0.25, 0.5 >, the 

distribution on arc (A,B) 2 is uniform on [0,3], and 
the distribution on arc (A,B) 3 is uniform on [2,3]. 
Note that the uniform distribution was used in many 
examples simply for ease of illustration. In general, 
the same reduction techniques apply for any discrete 
distribution. 


a (A Bji - <2,3,5> qjI — <.25,.25,.5> 



The difference between calculating the shortest 
and longest path lengths through a network mani- 
fests itself in reducing the parallel construct. For 
any given realization of a stochastic network that 
contains this construct, the shortest or longest path 
contains at most one of the parallel arcs shown in 
figure 5. When finding the shortest (longest) path 
through such a deterministic realization, only an arc 
having the minimum (maximum) length among these 
parallel arcs is considered for the shortest (longest) 
path. It follows that the set of parallel arcs can be 
replaced by a single arc (A,B), whose arc lengths rep- 
resent the minimum (maximum) lengths from all pos- 
sible combinations of arc lengths of the parallel arcs. 
One possible combination of arc lengths in the exam- 
pie is L{ A, B) 1 = 5, L( A, B) 2 = 2, and L{ A, B) 3 = 3. 
When reducing this construct with respect to finding 
the shortest path length, min(5,2,3) = 2 is included 
in G(a,B)- For this particular combination, 

P[P(A, B) 1 = 5]P[L(A,B) 2 = 2]P[L(A,B) 3 = 3] 

- 0.0625 

which contributes toward the overall probability 
P[i(A, B) = 2], Specifically, when considering the 
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shortest path length for this example, 


P[L(A, B) = 2] = £ P[L( A, B) 1 = 2}P[L(A, B) 2 = j]P[L( A, B) 3 = fc] 
j, k>2 


+ ^ P[L(A, B) 1 = j]P[L( A, B) 2 = 2]P[L(A, B) 3 = fc] 

J>2 

fc> 2 


+ ^ P[i(A,B) 1 =i]P[L(A,B) 2 = /c]P[L(A,B) 3 = 2] 
j, fc>2 

= 0.125 + 0.1875 + 0.09375 = 0.40625 

again using the statistical independence of arc lengths. When considering the longest path length through this 
example, 

P[L(A, B) = 2] = 52 WA, B) 1 = 2]P[L(A, B) 2 = j]P[L(A, B) 3 = k] 
j, k < 2 


+ 52 P[L( A, B) 1 = j}P[L( A, B) 2 = 2]P[L(A, B) 3 = k] 

3< 2 
k<2 


+ 52 P[A(A, B) 1 = j]P[L(A, B) 2 = k]P[L(A, B) 3 = 2] 
j , A;<2 


= 0.09375 + 0 + 0 = 0.09375 


More generally, any length in the set G(a,B) that 
arises from several combinations is represented only 
once in G(a,B) an< ^ corresponding probability is 
the sum of the appropriate individual probabilities. 
By applying this strategy, the parallel construct in 
figure 5 can be replaced by one of the single arcs 
shown in figure 6, depending on whether the shortest 
or the longest path length is of interest. 


reduced two-node network. For example, consider 
the network in figure 7, where the distribution on 
each arc is uniform on [1,2]; that is, = <1,2> 

and P(ij) = <0.5, 0.5 > for each (i,j) in the graph. 
By applying only series and parallel reductions, this 
network can be reduced to the equivalent shortest 
path network shown in figure 8. 


The series and parallel reductions are well-known 
and have applications in many optimization prob- 
lems. Martin (ref. 9) was the first to apply these 
reductions to the stochastic shortest path problem. 
By using these two reduction steps, more complex 
networks can sometimes be reduced to a single arc 
between the source and terminal nodes. In fact, 
a series-parallel network is defined to be any net- 
work that can be simplified via series and parallel 
reductions to an equivalent two-node network, that 
is, equivalent in the sense that the distribution of 
the shortest (or longest) path length for the origi- 
nal network is exactly the same distribution for the 


Distribution of 
shortest path length 


0 - 


a (A,B) = <0»1|2,3,5> 


-KB 


p (A B) = <.25 f . 25,. 40625,. 09375, 0> 


Distribution of 
longest path length 

3/* d\ = <0,1 ,2,3,5> 

0 — Kb) 

P(A B) = <0, 0,. 09375, .40625, .5> 
Figure 6. Reduced parallel constructs. 
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a/c t\ = <3,4,5, 6> 

© (s ' t) — ►© 

P( S , t ) = <-31 6406, .529297, .150879,. 00341 8> 

Figure 8. Reduced series-parallel network. 

The last fundamental reduction, called a figure- 
eight reduction, while less well-known than the se- 
ries and parallel reductions, incorporates the same 
approach to simplifying a network structure into an 
equivalent, less complex structure. The figure-eight 
construct, shown in figure 9, is found only in cyclic 
networks. In this construct, the center node B in 
the structure must have exactly indegree(B) = 2 and 
outdegree(B) = 2 and must be connected to exactly 
two other nodes, labeled A and C in figure 9. Specif- 
ically, B has one incoming arc from node A and one 
from C; both node A and node C have one incoming 
arc from B. 



the shortest path, the only remaining paths of inter- 
est are (A,B)^(B,C) and (C,B)—(B,A). The figure- 
eight construct is then simplified by replacing arcs 
(A,B) and (B,C) with a single arc (A,C), and by re- 
placing arcs (C,B) and (B,A) with a single arc (C,A). 
The lengths and probabilities of the new arcs are de- 
termined in the same manner as with series reduc- 
tions. For example, the set of lengths for (A,C) is the 
set of all possible additive combinations of lengths 
from (A,B) and (B,C). The simplified figure-eight as- 
sumes the form shown in figure 10. 



Figure 10. Reduced figure-eight construct. 

The advantage of using these fundamental reduc- 
tions is that they yield a single equivalent network 
that has a smaller number of arcs and possibly fewer 
nodes. Hence, determining the distribution of the 
shortest path length for this reduced network is eas- 
ier. In fact, these fundamental reductions have been 
previously applied in calculating network reliability 
(ref. 10). However, many complex network archi- 
tectures cannot be simplified via series and parallel 
reductions. The contribution made to the stochastic 
shortest path problem through this investigation is 
the development of a new reduction technique, called 
conditional factoring, that, like the fundamental re- 
ductions, is based on the configuration of the nodes 
and arcs in the network. Conditional factoring of- 
fers additional possibilities for simplifying a given 
network. 

Conditioned Reductions 


Figure 9. Figure-eight construct. 

Since the figure-eight construct occurs only in 
cyclic networks, only the shortest path from the 
source node to the terminal node is of interest. In 
this construct, four possible subpaths exist through 
node B: (A,B)^(B,A), (B,CMC,B), (A,B)-> 
(B,C), and (C,B)— ►(B,A). Any source-to-terminal 
(s-t) path through the network that contains the sub- 
path (A,B)-+(B,A) must be at least as long as the 
shortest path through the network. Since this sub- 
path does not affect the length of the shortest path, 
this subpath is irrelevant to the distribution of the 
shortest path length and can be ignored. Similarly, 
the subpath (B,C)— >(C,B) can safely be ignored. By 
eliminating these two paths from the set of possi- 
ble paths through node B that could be included in 


For a stochastic network whose architecture can- 
not be simplified by applying series and parallel 
reductions, most exact methods for finding the dis- 
tribution of the shortest path length depend on com- 
plete state-space enumeration (ref. 2) or enumeration 
of all cutsets or possible paths (refs. 7, 11, and 12). 
In simplifying a construct by complete enumeration, 
all possible combinations of arc lengths for all arcs in 
that construct are considered. For example, given 
the network construct in figure 11, if \(i,j)\ — 3 
for each of the 7 arcs in that construct, 3 7 = 2187 
states of that construct would be considered in solv- 
ing the network. If | (i, j) | = 10 for each arc in that 
construct, 10 7 = 10 million states of that network 
are considered. In general, complete enumeration 
is computationally infeasible except for very small 
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networks where the arcs can assume only a minimal 
number of lengths. 



Figure 11. Partial network structure. 

In this investigation, some constructs that are 
prevalent in many network architectures were iden- 
tified that can be simplified without complete enu- 
meration. The technique used to simplify these con- 
structs, called conditional factoring , is based on the 
idea that it is sufficient to consider only a specific sub- 
set of arcs within certain network constructs in order 
to generate an equivalent, simplified representation. 
Thus, the conditional reductions can significantly de- 
crease the computational requirements necessary to 
solve many networks and, hence, allow the exact so- 
lution of a larger class of networks. For example, if 
|(i,j)| = 3 for each arc in figure 11, only 9 states of 
that construct are needed to solve the network using 
conditional factoring; and, if |(z, j)| = 10, only 100 
states are necessary. In general, conditional factor- 
ing can be applied to any node within the network 
except the source or terminal nodes. However, con- 
ditional factoring reduces the computational effort 
only when it is applied to the four special constructs 
described in this section. For all other constructs, 
conditional factoring is believed to be equivalent to 
complete enumeration. 

The concept of conditional factoring is similar to 
the factoring theorem, which is used for computing 
network reliability (ref. 13), except that factoring 
is based on a specific node in the graph instead of 
an arc. The conditioned reductions are performed 
by identifying a specific construct centered around a 
node, called the central node, in the network. The 
structure of this construct is simplified by removing 
the central node and completing all possible paths 
through that central point of that construct. In con- 
trast to the fundamental reductions, conditional fac- 
toring generates a set of independent subnetworks to 
be solved, where each subnetwork incorporates the 
new, reduced structure into the overall network ar- 
chitecture. These subnetworks are generally smaller 
in size (i.e., have fewer nodes) and involve less compu- 
tational effort to solve. The lengths of the new arcs in 
the subnetworks are defined by conditioning on the 
arcs that are used more than once in defining new 
arcs in the reduced construct. These arcs that are 


used more than once in forming the new arcs of the 
reduced construct are called factoring arcs. A unique 
subnetwork is defined for each combination of lengths 
among the factoring arcs, and this eliminates any de- 
pendencies among the subnetworks. In this section, 
the concept of conditional factoring is illustrated by 
applying it to four special network constructs. 

The first construct to be introduced is the fan 
construct. A deterministic fan construct with central 
node B is shown in figure 12. Note that in any 
deterministic construct, the length of each arc is 
known with certainty; that is, the arc values are 
not random variables. In general, a fan construct 
with central node B has the following characteristics: 
(1) B is any node in the network except the source 
or terminal node, and (2) either indegree(B) = 1 or 
outdegree(B) = 1 (if both are 1, then B is the central 
node in a series structure). If indegree(B) = 1, 
the incoming arc to node B is called the factoring 
arc (since it is the only arc in the construct used 
more than one time in reducing the construct) and 
is denoted /(B). Similarly, if outdegree(B) = 1, the 
outgoing arc from node B is called the factoring arc 
and is denoted /(B). In figure 12, the factoring arc 
is /(B) = (A,B). 



Figure 12. Deterministic fan construct. 

In a deterministic construct, each arc (i,j) has 
a constant length, say k(i,j). Thus, within that 
construct, |(i, j)\ = 1 and P[L(i,j) = fc(i,j)] = 1. 
If two arcs, (A,B) and (B,C), with constant lengths 
are combined, the length of the resulting arc is the 
sum of the two arc lengths. Thus, the new length 
L(A, C) = k( A, B) + A;(B, C) is also a constant and 
has corresponding probability of 1. 

A fan construct is graphically simplified by elimi- 
nating the central node and completing each possible 
path (relative to the graph) through that central con- 
nection point. In a deterministic fan construct, the 
length of the factoring arc is added to the lengths 
of each of the other arcs connected to the central 
node. For example, arc (A,C) in figure 13 is formed 
from arcs (A,B) and (B,C), arc (A,D) is formed from 
arcs (A,B) and (B,D), and arc (A,E) is formed from 
arcs (A,B) and (B,E). The length of each new arc 
is accordingly the sum of lengths of its component 
arcs; and, since each arc has a constant length, the 
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probability associated with each new arc is 1. For the 
network in figure 12, this reduction process yields the 
network in figure 13. 



P[L(A,C)=6] = 1.0 


P[L(A,D)=7] = 1.0 


P[L(A,E)=10] = 1.0 


Figure 13. Reduced deterministic fan construct. 

The concept of conditional factoring applies when 
the length of each arc is assumed to be a random vari- 
able. To illustrate, consider the same graph shown in 
figure 12 but allow each arc to be a random variable. 
Figure 14 shows this construct where each arc can 
take on values in the n-tuple <1, 2> with correspond- 
ing probabilities < 0.4, 0.6 >. Such a stochastic fan 
construct centered at node B is decomposed based on 
the distribution of lengths on the factoring arc /(B). 
In figure 14, /(B) = (A,B) just as in figure 12. 


a (A,B) = a (B,C) = a (B,D) = a (B,E) = <1 ’ 2> 
P(A,B) = P(B,C) = P(B,D) = P(B,E) = <4 ” 6> 



Figure 14. Stochastic fan construct. 

Since the length of /(B) is now a random variable, 
a distinct subnetwork is generated for each length in 
a /(B) ■ The number of subnetworks created by fac- 
toring on node B is denoted w(B). For the stochastic 
fan construct, the number of subnetworks generated 
is equal to the number of possible lengths for the 
single factoring arc; that is, u;(B) = |/(B)|. To gen- 
erate each subnetwork Sj, for j = 1, 2, ..., w(B), 
the following steps are taken: (1) let /(B) take on 
a constant length Cj, where cj is the jth length in 
°/(B)i ( 2 ) eliminate node B from the original con- 
struct; (3) complete all possible paths through the 
central point; and (4) define the new arc lengths and 
the associated probabilities. 

In general, the arc lengths in the reduced con- 
struct are determined by adding the lengths of their 
component arcs. In the strictly deterministic con- 
struct, this is simple since only constant length arcs 
are combined. However, in reducing a stochastic 
construct, constant length axes (resulting from the 


factoring arc) are combined with arcs with random 
lengths. When an arc with a constant length is com- 
bined with an arc whose length is governed by a dis- 
crete distribution D, the length of the resulting arc 
is a random variable governed by D. 

The fan construct shown in figure 14 can be de- 
composed into an equivalent set of two independent 
subnetworks, shown in figure 15, since |/(B)| = 2. 
The first subnetwork S\ is constructed by assuming 
L(A, B) = 1, and the second subnetwork S 2 is con- 
structed for L(A, B) = 2. The probability of obtain- 
ing subnetwork 1 is P(Sj) = P[L(A, B) = 1] = 0.4, 
and, similarly, the probability of obtaining subnet- 
work 2 is P{S 2 ) = P[L{ A, B) = 2] = 0.6. 

The next three reductions focus on constructs 
with a central node B where either indegree(B) = 2 
or outdegree(B) = 2. Additionally, each of these 
constructs contains at least one simple cycle of size 
two (a cycle involving only two nodes) involving the 
central node B. 

The first of these constructs, shown in fig- 
ure 16, is called the loop construct with central 
node B. In general, a loop construct with central 
node B has the following characteristics: (1) inde- 
gree(B) = 2, (2) outdegree(B) = 2, and (3) B and 
some other node D form the only cycle of the con- 
struct, which is a simple cycle. The possible paths 
through this structure that are considered for the 
shortest path are then (A,B)— >(B,D), (A,B)— >(B,C), 
and (D,B)— >(B,C). The only other path through 
the structure, (A,B)— >(B,D)— »(D,B)— >(B,C) is never 
considered as part of a shortest path since a subset of 
that path, (A,B)— >(B,C), is always at least as short. 
As a result, we can replace each of these three paths 
with an equivalent arc, as shown by the reduced con- 
struct shown in figure 17. 

Since the arcs (A,B) and (B,C) are used more 
than once in defining the new arcs in the reduced 
construct, the lengths of both of these arcs must 
be considered when generating subnetworks. Hence, 
arcs (A,B) and (B,C) are factoring arcs for this struc- 
ture and are denoted /i(B) and / 2 (B), respectively. 
In general, the number of subnetworks produced by 
any conditioned reduction is the product of the dis- 
tribution sizes for each factoring arc. For a loop con- 
struct, a distinct subnetwork is produced for each 
combination of lengths for the two factoring arcs, so 
W (B) = |/i(B)||/ 2 (B)|. To generate the subnetwork 
S k , where k = 1, 2, ...,w(B), the following steps are 
taken: (1) let the factoring arcs /i(B) and / 2 (B) 
take on one of the combinations of constant values 
from Ofr(B) and ay 2 (B); (2) eliminate node B from 
the original structure; (3) complete all possible paths 
through the central point; and (4) appropriately de- 
fine the new arc lengths and their probabilities. 
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a (A,B) = a (B,C) = a (B,D) = a (B,E) = ^ - 2> 
P(A,B) = P (B,C) = P(B,D) = P(B,E) = < - 4 - 6> 



Figure 15. Subnetworks for the stochastic fan construct. 




Figure 17. Reduced loop construct. 

To illustrate, suppose that in the loop construct 
of figure 16 each arc can take on lengths in the 
n-tuple <1,2> with equal probability. From the first 
reduction step, four distinct subnetworks are gener- 
ated where L{ A,B) = 1 and L(B, C) = 1 are in Si, 
L(A,B) = 1 and L( B,C) = 2 are in S 2 , L(A,B) = 2 
and L(B,C) = 1 are in S3, and L(A,B) = 2 and 


L( B,C) = 2 are in S4. These four subnetworks 
are shown in figure 18. To simplify the notation in 
the following figures, the lengths and probabilities 
are displayed in an abbreviated form. In this nota- 
tion, the lengths and corresponding probabilities for 
each arc appear in the format: <length, probability: 
length, probability: ... : length, probability>. 

The third construct is a general loop construct 
and is shown in figure 19. The general loop construct 
with central node B has the following characteristics: 
(1) either indegree(B) = 2 and outdegree(B) > 2 or 
outdegree(B) = 2 and indegree(B) > 2; and (2) B 
together with some node A forms the only cycle of 
the construct, which again is a simple cycle. When 
the central node is removed from a general loop con- 
struct, the number of arcs in each resulting sub- 
network is at least as great as the number of arcs 
in the original network. For a general loop con- 
struct with central node B, if indegree(B) + out- 
degree^) = 5, the number of arcs in the reduced 
construct is the same as in the original construct. If 
indegree(B) + outdegree(B) > 5, the reduced con- 
struct contains more arcs than the original. 
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Figure 18. Subnetworks for the reduced loop construct. 

The reduction of a general loop construct is illus- 
trated in figure 20. In this case, the reduced con- 
struct has one more arc than the original construct. 
Note that in reducing the general loop construct each 
arc in the original structure is used more than once 
except for arc (B,A). As a result of this commonality, 
each arc in the original general loop construct must 
be a factoring arc except for arc (B,A). The number 
of factoring arcs 7 is thus 7 = indegree(B) + out- 
degree(B) - 1 and the number of subnetworks gen- 
erated is cj(B) = |/l(B)||/2(B)| . . . |/ 7 (B)|. As in the 
fan and loop constructs, the basic steps to gener- 
ate subnetwork 5^, where k = 1, 2, ... o;(B), are 
as follows: (1) let the factoring arcs fi{ B), where 
i = 1, 2 ,..., 7, take on one of the combinations of 
constant lengths in ay.(B); (2) eliminate node B from 
the original structure; (3) complete all possible paths 
through the central point; and (4) appropriately de- 
fine the new arc lengths and probabilities. 



indegree(B) = 2 



Figure 19. General loop constructs. 


Suppose each arc in the general loop structure of 
figure 20 can assume an integer length in the n-tuple 
<1,2> with equal probability. Then, there would 
be 7 = 5 factoring arcs, and u;(B) = 2 5 = 32 sub- 
networks generated in reducing that structure, com- 
pared with 2 6 = 64 networks that would be gener- 
ated using complete enumeration. Two of the subnet- 
works are shown in figure 21. In these subnetworks, 
only arc (C,A) has a nonconstant distribution. This 
is because (C,A) is formed by combining arcs (B,A) 
and (C,B), and arc (B,A) is the only arc in the struc- 
ture that is not a factoring arc of node B. 

The last conditioned reduction construct, shown 
in figure 22, is called a double loop construct. Like 
the general loop construct, either the indegree or 
the outdegree of the central node must be 2. The 
double loop construct, though, must contain exactly 
two simple cycles connected to the central node. 
In general, a double loop construct with central 
node B has the following characteristics: (1) either 
indegree(B) = 2 and outdegree(B) > 2, or out- 
degree^) = 2 and indegree(B) > 2; and (2) B forms 
exactly two simple cycles: one with some node A and 
one with some other node C; see figure 22. 
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When the double loop construct is reduced, two 
of the arcs in the original network are used only one 
time, such as arcs (B,A) and (B,C) for the construct 
with indegree(B) = 2 and arcs (A,B) and (C,B) for 
the construct with outdegree(B) = 2 in figure 22. As 
with all conditioned reductions, each arc that is used 
more than once in creating the reduced construct is 
a factoring arc. Hence, for any double loop construct 
with central node B there are 7 = indegree(B) + out- 
degree(B) - 2 factoring arcs, and the total number 
of subnetworks is u;(B) = |/i(B)||/ 2(B)| . . . |/ 7 (B)|. 
The steps used to generate each subnetwork of a 
double loop construct are the same as those used 
for the other conditioned reductions. An example of 
the graphical reduction of a double loop construct 
is shown in figure 23. As with the general loop 
construct, removing the central node of a double 
loop construct can produce more arcs in each result- 
ing subnetwork than are present in the original con- 
struct. When indegree(B) + outdegree(B) > 6, more 
arcs are present in the reduced construct than in the 
original. Some of these additional arcs, though, may 
lead to further parallel reductions. 




Figure 20. Reduction of general loop construct. 




Figure 21. Two subnetworks for the reduced general loop 
construct. 

Conditioned Reductions and Shortest 
Path Calculations 

The rules for simplifying a special construct 
present within a stochastic network via conditional 
factoring follow a general scheme. Namely, once one 
of the constructs discussed in the preceding section 
has been identified, the following steps are taken to 
generate the desired subnetworks: 

1. Identify all factoring arcs for that construct. 

2. List all possible combinations of lengths for 
the factoring arcs. Each combination defines 
a distinct subnetwork. 

3. For each combination of lengths, compute the 
product of the probabilities associated with 
each length in the combination. Because of 
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indegree(B) = 2 



outdegree(B) = 2 

Figure 22. Double loop constructs. 



Figure 23. Reduction of a double loop construct. 


the statistical independence assumption, this 
product is the probability of that subnetwork 
occurring. 

4. Remove the central node from the construct 
and graphically complete all possible paths 
through the construct. (Note that any loop 
arc (i,z) formed can be eliminated from that 
subnetwork structure.) 

5. Add the lengths of the component arcs to get 
the lengths of the new arcs in the subnetworks. 

6. Define the probabilities for the new arcs. An 
arc with a constant length has probability 1. 
When an arc with a constant length is com- 
bined with an arc whose length is governed 
by a discrete distribution D, the length of the 
resulting arc is a random variable governed 
by D. 

By following these steps, a network can be decom- 
posed into a set of subnetworks that are more easily 
solved. To demonstrate how these steps are used to 
calculate the exact distribution of the shortest path 
length, consider the network shown in figure 24; no- 
tice that this network is not a series-parallel network. 
Suppose each arc in figure 24 can take on lengths in 
the n-tuple <1,3> with equal probability. 



Figure 24. Network with loop construct. 


In this network, there is a loop construct with 
central node B and a loop construct with central 
node C. For this example, the choice of central node 
is not important, so node B is arbitrarily selected. 
The heuristic rules used to govern the choice of 
which node to remove are discussed in a later section. 
According to steps 1 and 2 of the decomposition 
steps, there will be four subnetworks generated from 
the four possible combinations of arc lengths from the 
factoring arcs (A,B) and (B,D). These combinations 
are L(A,B) = 1 and L(B,D) = 1 in Si, L(A, B) = 1 
and L( B, D) = 3 in S2, L( A, B) = 3 and L( B, D) = 1 
in £3, and L(A, B) = 3 and L(B,D) = 3 in S4. 
The four subnetworks are shown in figure 25. From 
step 3 of the decomposition steps, the probability 
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associated with each subnetwork is the product of the 
probabilities corresponding to the constant lengths 
used to define that subnetwork. In figure 25, 

P(Si) = P[L(A,B) = 1]P[L(B,D) = 1] = 0.25 
P(S 2 ) = P[L( A,B) = 1]P[L(B,D) = 3] = 0.25 
P(S Z ) = P[L{ A,B) = 3]P[L(B,D) = 1] = 0.25 
P(S\) = P[L(A,B) = 3]P[L(B,D) = 3] = 0.25 


< 2 , 1 , 0 > 

'^2,.5:4,.5> 

\ \<2,.5:4,.5^ 

<1,.5:3,.5> V<1,.5:3,.5> 

Si 



<4,1 .0> 

\ \<4,.5:6,.5>/ 

<1,.5:3,.5> V<i,.5:3,.5> 

s 2 



<4,1 ,0> 



< 6 , 1 . 0 > 



<1 ,.5:3,.5> ^<i > .5:3,.5> 

s 4 


Figure 25. Subnetwork structures for the reduced network. 


Notice that each subnetwork in figure 25 can now 
be simplified to an equivalent two-node network, with 
respect to finding the distribution of the shortest 
path length, via series and parallel reductions. These 
simplified networks are shown in figure 26. 


^ ^ <2,1 .0> ^ 

s, ® *■<§) 

<2,.25:3,.125:4,.625> 

s 2 ® — » © 

<2,.25:3,.125:4,.625> 

S3 ® — * ■© 


_ „<2,.25:4,.5:6,.25>^-. i 

s 4 ® 


Figure 26. Reduced subnetworks. 


The distribution of the shortest path length 
through the original network in figure 24 is calculated 
as shown in table 1. The distribution for the origi- 
nal network is found by combining the shortest path 
distributions in each subnetwork Si using as weights 
the appropriate P(Si)‘, that is, for each subnetwork, 
the probability of each path length is multiplied by 
the probability of that subnetwork occurring and the 
resulting probabilities are summed across all subnet- 
works. Note that the columns for the distributions in 
table 1 all sum to 1.0 since the shortest path length 
distribution for each individual subnetwork must sum 
to 1.0. 


Table 1. Distribution of Shortest Path Length in the Network With a Loop Construct 


Length of 
shortest path 

Probability in final subnetworks 

Probability for 
original network 

Si 

P(Si) = 0.25 

s 2 

P{S 2 ) = 0.25 

S3 

P(s 3 ) = 0.25 

Sa 

P(S 4 ) = 0.25 

2 

1.0 

■ 




3 

0 





4 

0 





6 

0 

0 

0 


.0625 
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Total Factoring 

Although the special constructs required for the 
basic and conditioned reductions are found in many 
network configurations, realistic network configura- 
tions often do not, on initial inspection, contain any 
of these constructs. In these cases, the factoring ap- 
proach and the steps for conditioned reductions can 
still be applied; however, all arcs that are incident 
with the node chosen for removal are now factoring 
arcs. That is, all possible combinations of arc lengths 
around a central node are completely enumerated. 
If total factoring is required for each node of the 
network except for the source and terminal nodes, 
this approach is equivalent to complete enumeration, 
which, as discussed earlier, is computationally infea- 
sible for all but the smallest networks. Fortunately, 
applying total factoring on a node in a complex net- 
work often yields subnetworks that can be simplified 
via the basic and conditioned reductions. 

For example, the network shown in figure 27 does 
not contain any of the special constructs discussed in 
the previous sections. 



Figure 27. Network with no special constructs. 

To demonstrate total factoring, let node B be the 
central node to be removed from the network. The 
subnetworks generated by removing node B each 
have the form shown in figure 28. 



Figure 28. Subnetwork structure after total factoring. 

After performing a parallel reduction on the arcs con- 
necting nodes A and C, we can now identify a gen- 
eral loop construct with central node D in this sub- 
network structure. When node D is removed, the 
new subnetwork shown in figure 29 results. In this 


new subnetwork, there is a simple loop construct with 
central node C, a general loop construct with central 
node F, and a double loop construct with central 
node E. Since there are fewer factoring arcs associ- 
ated with the simple loop construct than with the 
general or double loop constructs, node C is chosen 
for factoring. After factoring on node C, the sub- 
networks whose structure is shown in figure 30 are 
generated. After performing all appropriate parallel 
reductions, these subnetworks have the same graphi- 
cal form as the example given in figure 24 and, thus, 
can be solved in the same manner. 



Figure 29. Second subnetwork structure after reduction of 
general loop construct. 



Figure 30. Third subnetwork structure after reduction of 
simple loop construct. 

Determining the Shortest Path Length 
Distribution 

In previous sections, the distribution of the short- 
est path length has been found for series-parallel net- 
works and for networks where only one level of sub- 
networks needs to be generated. For many networks, 
such as the one given in figure 27, reductions must be 
repeatedly applied before the original network is sim- 
plified to an equivalent two-node network. That is, 
each generated subnetwork can itself be reduced by 
further application of conditioned and fundamental 
reductions. 

In general, the final subnetworks that simplify 
to two-node structures (weighted by their associated 
probabilities of occurrence) are the only contributors 
to the overall distribution of the shortest path length 
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in the network. Hence, for each possible subnetwork, 
it is important to keep track of the probability of 
that subnetwork occurring. Suppose a conditioned 
reduction is applied to a network G, and several 
subnetworks, denoted 5^, are generated. Then, by 
the statistical independence assumption, the proba- 
bility of each subnetwork is calculated according 
to step 3 in the general reduction procedure. Now 
suppose that another conditioned reduction is used 
to simplify subnetwork 5^ and additional j subnet- 
works, denoted are generated. The current sta- 
tus of the computations can be represented by the 
tree depicted in figure 31. The shaded leaves of this 
subnetwork tree represent the currently unresolved 
subnetworks. 



Figure 31. Subnetwork tree. 

The probability of the j th subnetwork of 5& occur- 
ring is P[S kJ ] = P[S kJ \S k ]P[S k ], where P[S kJ \S k ] 
is calculated according to step 3 in the decompo- 
sition procedure. This process of generating new 
subnetworks from existing subnetworks continues un- 
til a generated subnetwork simplifies to a two-node 
structure. 

Note that all subnetworks at level 2 of the sub- 
network tree have the same graphical structure. The 
only differences among these subnetworks are the 
lengths and probabilities occurring on the arcs. More 
generally, all subnetworks on the same level of the 
subnetwork tree have the same graphical structure. 
So if a conditioned reduction is required to sim- 
plify Sfc, the same conditioned reduction applies to 
*?l, 52, 5fc_ i- Also, if some S^j in figure 31 simpli- 
fies to a two-node network without additional condi- 
tioned reductions, each subnetwork on level 3 also 
simplifies to a two-node network without generat- 
ing further subnetworks. These characteristics of the 
subnetwork tree suggest that the reduction algorithm 
may be amenable to parallel implementation. 

To calculate the distribution of the shortest path 
length in the original network, only the subnetworks 
at the final level of the subnetwork tree are consid- 
ered; that is, those leaf subnetworks that simplify to 
two-node networks without further conditioned re- 
ductions. For each leaf subnetwork in the final level, 
the probability of each arc length appearing in the 
associated distribution is multiplied by the proba- 


bility of that subnetwork occurring. Then, for ev- 
ery possible length, the contributing probabilities are 
summed across all leaf subnetworks. This yields the 
distribution of the shortest path length in the original 
network. 

This process of calculating the distribution of 
the shortest path length has been implemented in 
a computer program. The program is written 
in FORTRAN 77 and is implemented on a Digi- 
tal Equipment Corporation VAXstation 3200 work- 
station. The algorithm for determining the distribu- 
tion of the shortest path length is given as follows: 

Input network data and store it 

Use all applicable basic reductions to simplify the 
network 

If the network has been completely reduced 
then 

Store the distribution information (note: 
this is the complete distribution of the 
shortest path length) 

else 

Identify a node to factor on 

Apply the appropriate reduction, generate 
the subnetworks, and place them in a 
stack 

While there are subnetworks in the stack 

Remove the subnetwork on the top of 
the stack 

Use basic reductions to simplify the 
subnetwork 

If the subnetwork has been completely 
reduced 

then 

Store the distribution informa- 
tion for that subnetwork 

else 

Identify a node within the sub- 
network to factor on 

Apply the appropriate reduc- 
tion, generate new subnetworks 
and place them on top 
of the stack 

end if condition 
end while condition 
end if condition 

Print the distribution of the shortest path length. 
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The algorithm used to find the distribution of the 
longest path length is identical except that maximum 
in place of minimum values are used in computing the 
parallel reductions. 

In theory, this algorithm is valid for an arbitrar- 
ily large network. However, the implementation of 
the algorithm is limited since the number of subnet- 
works that must be generated and stored to solve a 
problem can overwhelm the computational resources 
available. In general, the size of the subnetwork tree 
depends on the complexity, number of nodes, and 
size of the distribution of lengths for each arc in the 
original network. If there are n nodes in the origi- 
nal network, there will be at most n — 1 levels in the 
subnetwork tree. If there are k subnetworks gener- 
ated during each simplification, then there would be 
a total of 1 + k + k 1 + ... + k n ~ 2 = ( k n ~ l - 1 )/(k — 1) 
subnetworks to solve. However, the maximum num- 
ber of subnetworks that ever need to be stored is 
(k— l)(n— 1) — (k— 2). These are worst-case estimates. 

In implementing the algorithm, an effort was 
made to minimize the size of the subnetwork tree by 
judiciously choosing the node for factoring. Specifi- 
cally, the factoring node i is chosen to fulfill the fol- 
lowing criteria: 

1. For all nodes j in the current subnetwork, 
min[indegree(i), outdegree(z)] < min[indegree 
(j), outdegree(j)]. 

2. Among all nodes k that satisfy condition 1, 

cj(z) < cj(fc). 

3. Among all nodes that satisfy conditions 1 and 
2, i is the first node that occurs in the graph. 

Empirical evidence has shown that in many cases, 
choosing the factoring node according to these crite- 
ria results in a smaller total number of subnetworks 
that need to be solved. However, these criteria are 
not optimal for all networks. Further research is nec- 
essary to determine whether globally optimal criteria 
exist for choosing the order of factoring. In practice, 
the number of subnetworks solved and the maximum 
number stored have been far fewer than those pre- 
dicted by the worst-case estimates. The results pre- 
sented in the next section demonstrate that a fairly 
modest computational effort is usually required. 

Applications of the Structural Factoring 
Approach 

In this section, the structural factoring approach 
is applied to several network problems taken from 
the literature. Two of the more interesting problems 
associated with network analysis are determining the 
distribution of the longest, or critical, path and the 


two-terminal reliability of a network. These two 
problems are actually special cases of the stochastic 
shortest path problem. 

Critical Path Analysis 

Although the analysis of the critical path is lim- 
ited to acyclic networks, critical path analysis has 
many applications to scheduling and performance 
problems such as those associated with communica- 
tions networks. To determine the distribution of the 
length of the critical path, the structural reductions 
are applied to an acyclic network just as they are 
in determining the shortest path except for decom- 
posing a parallel construct. Recall that in reduc- 
ing a parallel construct, the maximum arc lengths 
are considered instead of the minimum arc lengths. 
For example, consider the network problem presented 
by Fulkerson (ref. 6). The network has a simple 
four-node architecture, shown in figure 32, where the 
length of each arc is uniformly distributed on the 
interval [0,2]. 



Figure 32. Fulkerson network (ref. 6). 

In this network, two fan constructs can be identi- 
fied, one centered at node B, since indegree(B) = 1, 
and the other centered at node C, since outde- 
gree(C) = 1. The distribution of the longest path 
length through the network has been determined us- 
ing the implementation of the algorithm and is given 
in table 2. From the distribution in table 2, the ex- 
pected length of the critical path is calculated to be 
3.32510288. In the Fulkerson paper, the lower bound 
calculated by Fulkerson’s method is given as 3.22. 
This problem required 3.5 seconds of central process- 
ing unit (cpu) time to solve for the exact distribution 
and four subnetworks were generated. 

Another renowned problem, shown in figure 33, is 
the crossing network, also referred to as the “wheat- 
stone bridge” network, analyzed by Kleindorfer 
(ref. 14). Since this network is acyclic, the distri- 
butions for the length of the shortest and longest 
paths can be calculated. In his paper, Kleindorfer 
gives bounds on the cumulative distribution for this 
network where each arc assumes a length in the set 
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<1,2,3,4,5> with equal probability. By apply- 
ing structural factoring to the fan structures in this 
network, the exact distribution of the shortest path 
length from node A to node F can be determined. Ta- 
ble 3 gives the the distribution of the shortest path 
for this network; table 4 shows the cumulative dis- 
tribution of the critical (longest) path through the 
same network as well as the Kleindorfer bounds. 



Figure 33. Kleindorfer crossing network (ref. 14). 


Table 2. Distribution of Longest Path Length 
for the Fulkerson Network 


Length of 
longest path 

Exact probability 

0 

1/243 

1 

11/243 

2 

49/243 

3 

74/243 

4 

72/243 

5 

27/243 

6 

9/243 


Table 3. Distribution of Shortest Path Length 
for the Crossing Network 


Length of 
shortest path 

Probability 

3 

0.03064064 

4 

.08365312 

5 

.14335488 

6 

.18986496 

7 

.20426496 

8 

.16326144 

9 

.10479360 

10 

.05362176 

11 

.02052864 

12 

.00505344 

13 

.00087552 

14 

.00008448 

15 

.00000256 


Table 4. Distribution of Critical Path Length for the 
Crossing Network Compared With the Kleindorfer 
Cumulative Bounds 


Length of 
critical path 

Kleindorfer 
lower bound 
(ref. 14) 

Cumulative 

probability 

Kleindorfer 
upper bound 
(ref. 14) 

3 

0 

0.00000256 

0.008 

4 

0 

.00008704 

.032 

5 

0 

.00096256 

.080 

6 

0.002 

.00601600 

.160 

7 

.014 

.02654464 

.280 

8 

.055 

.08016640 

.424 

9 

.149 

.18496000 

.576 

10 

.312 

.34822144 

.720 

11 

.528 

.55248640 

.840 

12 

.731 

.74235136 

.920 

13 

.882 

.88570624 

.968 

14 

.969 

.96935936 

.992 

15 

1.000 

1.00000000 

1.000 


For critical path problems, the expected project 
duration time, which is the expected length of the 
critical path, is often of interest. The expected 
project duration was calculated from the exact dis- 
tribution of the critical path length to be 11.203136. 
Kleindorfer gives a lower bound (LB) of 9.000 and 
an upper bound (UB) of 11.358 for this problem. 
Shogan (ref. 15) also gives bounds on the expected 
project duration for this problem: LB = 10.6 and 
UB = 11.358. Only 3.7 seconds of cpu time were 
required to solve this problem exactly. 

Since each arc in this network has five possible 
lengths, the complete enumeration approach would 
consider 5 8 = 390625 possible states of the network 
to determine exactly the distribution of the shortest 
path length. For each possible state of the network, 
the shortest (or longest) path length must be identi- 
fied. In contrast, to solve this problem, the structural 
factoring approach generated only 31 subnetworks (of 
which 30 were less complex than the original net- 
work). A computer program that implements the 
complete enumeration approach took 57.21 seconds 
of cpu time to execute, and the resulting distribution 
agreed with the one given in table 3. The conditional 
factoring approach is clearly a drastic improvement 
over complete enumeration for determining an exact 
solution. 

Two-Terminal Reliability 

The translation of the stochastic shortest path 
problem to the two-terminal reliability problem is not 
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as obvious. The two-terminal reliability problem is 
defined relative to a network G in which each arc 
(z, j) has a given probability p(z,j) of functioning 
independent of all other arcs. The two-terminal 
reliability of the network is equal to the probability 
that there exists at least one path in G from a 
specified node s to a specified node t along which 
all arcs are functioning. 

Now consider a stochastic network G f involving 
the same nodes and same arcs as G. Each arc (z, j) 
in G f can assume only the lengths 0 and 1, with 
probabilities p(i, j) and 1 — p(z, j), respectively. If 
there is a shortest s-t (source to terminal node) path 
through G ' that has length 0, then there must be 
some path in G ' on which all arc lengths are 0. 
This directly corresponds to a path in G composed 
of all functioning arcs. Similarly, if the network 
is functioning, then there must be a path through 
the G f network that has length 0, and this will 
be a shortest s-t path. Thus, the probability that 
the shortest s-t path through G f has length zero is 
precisely the two-terminal reliability of the original 
network G. There are no architectural restrictions 
on the network when computing the two-terminal 
reliability except that there be a source and terminal 
node. Thus, the exact two-terminal reliability can 
be determined for cyclic as well as acyclic networks. 
Network reliability problems are known however to 
be #P-complete (ref. 1) where # P is a class of 
counting problems analogous to NP. Hence, #P- 
complete problems are at least as difficult to solve, if 
not more difficult than NP-complete problems such 
as the traveling salesman problem (ref. 1). 

The next example is a cyclic network, shown in 
figure 34, presented by Shogan (ref. 16). Two loop 
constructs, at nodes B and D, can be identified in 
the original network structure. 



For this example, the n-tuple of lengths for each 
arc is < 0, 2, 4 > with corresponding probabilities 
<0.4, 0.2, 0.4>. The exact distribution of the shortest 
path length through this network is given in table 5. 


Table 5. Distribution of Shortest 
Path Length for the Shogan Network 


Length of 
shortest path 

Probability 

0 

0.3047540654 

2 

.2810753843 

4 

.2784418939 

6 

.0986955899 

8 

.0339240387 

10 

.0028468838 

12 

.0002621440 


To solve this network for the distribution of the 
shortest path length, 1528 subnetworks were gener- 
ated, but the maximum number of networks stored 
by the computer at any given time was 27. Although 
a large number of subnetworks were generated, this 
problem took only 8.19 seconds of cpu time to solve. 
By contrast, it took 6048.89 seconds of cpu time 
to solve the same problem by completely enumer- 
ating all 2 15 = 32 768 possible combinations of path 
lengths. 

This problem can also be considered as a two- 
terminal reliability problem where the probability 
that each arc is operational is 0.4. The two-terminal 
reliability of the network, as shown in table 5, is 
0.3047540654. Shogan gives the following bounds on 
the two-terminal reliability: LB = 0.1971 and UB = 
0.3527. 

The last example, shown in figure 35, is a large 
network with 20 nodes and 38 arcs that was also pre- 
sented by Kleindorfer (ref. 14). The two-terminal 
reliability of this network, in which each arc is op- 
erational with probability 0.9, can be determined by 
letting each arc take on lengths < 0, 1 > with cor- 
responding probabilities <0.9, 0.1>. Using the pro- 
gram with these parameters, the s-t reliability of this 
network is found to be 0.98612801. The subnetwork 
tree generated for this problem had 15 levels. A max- 
imum of 15 subnetworks were stored at any given 
time, and 32 767 total subnetworks were generated. 
The program required 49.21 seconds of cpu time to 
solve this problem. For this problem, complete enu- 
meration would not be feasible since 2 38 or approxi- 
mately 2.75 x 10 11 states of the network would be in- 
dividually examined to find the shortest path length 
through each. Based on the time to solve a portion 
of this network by complete enumeration, approxi- 
mately 280 000 000 seconds of cpu time (or approxi- 
mately 9 years) would be required to solve the large 
Kleindorfer network. 
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Comparison Between Structural Factoring 
and Complete Enumeration 

The effort required to solve the examples pre- 
sented in the previous sections with structural fac- 
toring is compared in table 6 with the effort re- 
quired with complete enumeration. As the size of 
the network in terms of number of nodes and arcs 
increases, the total number of subnetworks and cpu 
time increases in an exponential manner as expected 
for an NP-complete problem. However, the struc- 
tural factoring approach is dramatically superior to 
the straightforward complete enumeration approach 


and allows exact solution to large network problems 
that were previously analyzed only by approximation 
techniques. 

Conclusions 

The structural factoring approach allows the ex- 
act distribution of the shortest path length through 
a stochastic network to be" determined. Applica- 
tion of structural factoring to fan, loop, general loop, 
and double loop constructs significantly increases the 
tract ability and applicability of the approach, since 
complete enumeration of all states is not necessary 
to simplify these constructs. Although there are net- 
work architectures that do not initially contain these 
special constructs, factoring can often be applied lo- 
cally to produce subnetworks that do contain such 
constructs. Hence, this approach extends the com- 
putational range, especially for cyclic networks, to 
larger and more complex networks than have previ- 
ously been solved with other exact methods. Struc- 
tural factoring can also be applied to find the distri- 
bution of the critical path length in acyclic networks 
and the two-terminal reliability for a network. 

Although the technique of structural factoring is 
theoretically unconstrained, the size of the network 
that can be solved with the algorithm is restricted 
by computational resources. In the worst case, the 
computational effort can grow exponentially with the 
size of the problem, which is indicative of the NP- 
complete characterization of the stochastic network 
problem. Computational experience with the struc- 
tural factoring algorithm has nonetheless been en- 
couraging. Moreover, the ability to solve certain 
nontrivial problems exactly can serve as a baseline 
for assessing the accuracy of various approximation 
schemes proposed to solve larger problems. 

NASA Langley Research Center 
Hampton, VA 23665-5225 
January 3, 1991 


Table 6. Comparison Between Structural Factoring and Complete Enumeration 


Network (n, m, k) a 

Structural factoring 

Complete enumeration 

No. of subnetworks 

cpu time, sec 

No. of states 

cpu time, sec 

Fulkerson (4,5,3) 

4 

3.5 

243 

1.03 

Kleindorfer crossing (6,8,5) 

31 

3.6 

390625 

57.21 

Shogan (7,15,3) 

1528 

8.19 

14 348 907 

6048.89 

Large Kleindorfer (20,38,2) 

8191 

49.21 

2.748779 x 10 u 

2.8 x 10 8 


n is the number of nodes in the network, m is the number of arcs, and k is the number of “lengths” for each arc. 
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